X and Y waves in the spatiotemporal Kerr dynamics of a self-guided light beam 
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The nonlinear stage of development of the spatiotemporal instability of the monochromatic Townes 
beam in a medium with self-focusing nonlinearity and normal dispersion is studied by analytical 
and numerical means. Small perturbations to the self-guided light beam are found to grow into two 
giant, splitting Y pulses featuring shock fronts on opposite sides. Each shocking pulse amplifies a 
co-propagating X wave, or dispersion- and diffraction-free linear wave mode of the medium, with 
super-broad spectrum. 



I. INTRODUCTION 

The loud consequences of the instability of the 
monochromatic Townes beam [l[ in self-focusing Kerr 
media under perturbations in space — catastrophic col- 
lapse or diffraction — , as well as the central role played 
by the Townes profile in the theory of self-similar col- 
lapse @, S 0] i have probably overshadowed its relevance 
in other contexts. The Townes beam suffers also from 
an instability under small spatiotemporal perturbations, 
which was predicted in [5| , but only studied in detail 
recently. [6| In its initial stages of growth, a linear 
instability analysis showed that the spatiotemporal in- 
stability is characterized by two weak, Y-shaped unsta- 
ble modes that can grow exponentially. These weak Y- 
shaped modes have been seen to account for the incipient 
spatiotemporal dynamics observed in a seemingly sta- 
tionary, self-guided light beam Q (e.g., a light filament), 
which includes the onset of spectral broadening and of 
conical emission, and the seed for temporal break-up [see 
for instance Ref. Q for these phenomena]. In numeri- 
cal simulations, weak Y-shaped wave modes have been 
seen to develop from spatiotemporal perturbations to a 
Townes beam. [6( In experiments, the single and dou- 
ble Y-shaped structures observed in the spatiotemporal 
spectra of light filaments immediately after collapse are 
compatible with the formation of the unstable modes of 
the newly formed self-trapped beam. [(J 

In this paper we take a step forward by studying the 
spatiotemporal instability of the Townes beam in its non- 
linear regime of development, where the amplitude of the 
perturbation becomes comparable to, or higher than the 
amplitude of the Townes beam, and the energy lost b 
the Townes beam is not negligible. Assuming from 
that the Y unstable modes have grown up to a small am- 
plitude, we show here that they develop into two giant, 
splitting (subluminal and superluminal) Y pulses that 
exhibit broadened spectra and shocks fronts on opposite 
sides. Moreover, the nonlinear interaction between each 
Y pulse and the Townes beam results in the exponen- 
tial amplification of two X waves with super-broad spec- 
tra that co-propagate with the shocking Y pulses. The 
Townes beam, the Y pulses and the X waves are seen to 
coexist during the propagation, interacting strongly in a 
highly phase-matched and group-matched configuration. 



The spatiotemporal dynamics of Y pulses and X waves 
are found to resemble very closely the observed dynamics 
of light filaments, from their formation up to their later 
stages of proapgation. @, fi ©, El EJ E d 

Thus, though mathematically complex, the theory on 
the spatiotemporal instability of a self-guided light beam 
in [6| (linear regime) and in the present paper (nonlinear 
regime) can explain from minimum ingredients (diffrac- 
tion, dispersion and the Kerr nonlinearity) the origin of 
the basic phenomena observed in the spatiotemporal dy- 
namics of a self-guided light beam, though their quanti- 
tative description would require, of course, the inclusion 
of all relevant higher-order effects, as higher-order dis- 
persion, plasma defocusing, nonlinear absorption, etc. 

After recalling briefly the linear instability analysis 
(Sec. HI|, we first study by analytical and numerical 
means the nonlinear regime of growth of a single Y per- 
turbation at a single frequency to the Townes beam (Sec. 
IIII|) . The results allows us to understand the most general 
case of the nonlinear development of the two unstable Y 
modes when excited at a continuous of frequencies (Sec. 

E2J. 



II. LINEAR SPATIOTEMPORAL INSTABILITY 
ANALYSIS 

We recall that the monochromatic Townes beam a = 
<Xa(p) exp(ia£), where ao(p) is the transversal Townes 
profile, and a ~ 0.2055 for a {p) = 1 [l|, is a transver- 
sally localized, self-guided solution to the cubic nonlinear 
Schrodinger equation (NLSE), which for a wave a(p, r, £) 
with revolution symmetry about the propagation direc- 
tion f, reads as 



1 d p (pd p a) 

2 P 



-d^a + i\a\ 2 a . 



(1) 



The term with the second derivative with respect to the 
local time r accounts for normal group velocity dispersion 
in the non-monochromatic case. Apart from the well- 
known instability of algebraical type under pure spatial 
perturbations, an instability under small spatiotemporal 
perturbations of exponential type has been described re- 
cently [(| . A couple of unstable modes of the form 



p K = e[u{p)e 



-%Q, s t-\-%k^ 



■v*(p)e 



iQ-sT — lK, 1 



(2) 



2 




FIG. 1: (a) (dashed curve) Gain — 2ki, and (solid curve) real 
part kr of the eigenvalue associated to the unstable modes 
p K and p- K * at each perturbation frequency ±Q S . (Dotted 
curve) Quotient kr/Q s , and (dot-dashed curve) the function 
+ a. (b) Modulus of the spatiotemporal spectrum of 
the unstable perturbation p K . Normalization is such that the 
power J °° \p K \ 2 pdp is the same for all Q s . 



and 



with kj = ImK < 0, for each temporal frequency shift 
±Q S with respect to the Townes beam frequency, can 
grow exponentially with gain —2kj. Figure [Ha) shows 
the gain — 2/tj, and kr = ReK > 0, that determines 
the wave number shift along the ^-direction of the un- 
stable modes with respect to the wave number of the 
Townes beam, for each frequency shift Q s of the pertur- 
bation. Figure [T^b) shows the spatiotemporal spectrum 
of the unstable wave mode p K , obtained by plotting the 
Hankel transforms u(Q) and v(Q) of u(p) and v(p) at 
their respective frequency shifts +fi s and — Q s [similar 
plot is obtained for p- K * by reflection of Fig. [TJb) about 
Q s = 0] . The structure of the spatiotemporal spectrum 
in the form of a vaulted Y-letter is due to the definite 
maxima at radial frequencies 



f for +fl s [i.e., for tt(Q)], 

H ~ \ v^|n s | for -n B [i.e., for v(Q)}. { > 

This structure points out a common origin for the ax- 
ial and conical radiation of new frequencies observed in 
self-guided light beams. Up-shifted new frequencies prop- 
agating collinearly with the Townes (Q ~ 0) beam grow 
jointly with new down-shifted, non-collinear frequencies 
(Q ~ \/2 | Sis |) as a consequence of the excitation of the 
unstable mode p K (and vice versa for p- K *). Y-shaped 
spectra have been shown in Q to develop spontaneously 
from perturbations in the Townes beam. In filamentation 
experiments, single and double Y-shaped spectra have 
been observed to emerge soon after collapse. [f| 



III. NONLINEAR SPATIOTEMPORAL 
INSTABILITY: SINGLE Y WAVE 
PERTURBATION 

A. Asymptotic analysis 

The validity of the above linear instability analysis is 
limited to initial stage of growth of the perturbations, 
where \p K \, \p- K * \ <C 1. Some properties of their further 
evolution can be inferred from the following analysis. As- 
sume, for simplicity, that only one of the two unstable Y 
modes, say p K , and at only one couple of frequencies ±Q S , 
is growing, as sketched in Fig. [2] (circles), starting with 
small amplitude e. For small enough £, the total field 
a = {a (p) + e[u(p)e- in ° T+iK t + v+ipy^'^^je 1 ^ is 
seen to satisfy the relation 

d^a — — kj e d e a — (Kn/fl s ) d T a , (5) 

where the reduced field a = aexp(— ia^) is introduced for 
convenience. Though this relation pertains to the small 
perturbation regime, or small enough £, it is not difficult 
to demonstrate that it continues to hold along the entire 
evolution up to £ — * oo (where the total field continues 
to depend parametrically on the starting amplitude e). 
As said, Eq. ||5j) is assumed to be satisfied at a certain 
small value of £, say £o, i-e. , 

<%a|? = [-Kied e a- (kr/O s ) S T a]|f . (6) 

From the NLSE the evolution equation for a, valid 
everywhere, is of the form d^a = C[a] + i\a\ 2 a, where £ 
is a linear differential operator. The evolution equation 
for a implies, in particular, that d^(d^a)\^ = d^[C[a] + 
i\a 2 \a]\ 5o . Making use of the linearity of £ and of relation 
© at £cb we find, after straightforward algebra, 

daid&ko = c%[-K/e<9 £ a- (kr/Q s ) <9 T a]| 5o • (7) 

The derivatives of the left hand and the right hand sides 
of Eq. ^ at £o are also equal. Analogously, starting 
from fl£(0£6)| fo = d£[£[a] + i\a 2 \a}\ io (n = 2,3.. .), we 
also find 

di( d ^)ko = dl l {- Kl ed e a-(K R /n s )d T a}\s . (8) 

Thus, the full nonlinear evolution governed by the NLSE 
([1]) implies that if relation ([6]) is valid at £o, it is also valid 
for all its derivatives at This assertion is equivalent 
to say that relation © is valid not only at £o but also at 
any propagation distance £. This last point can be seen 
in more formal terms by expressing d^a at any distance 
£ as the power series 

^ B = E 3 Wc a )l& (£-&)"• (9) 

Using Eq. ([5]) for the successive derivatives at £o in Eq. 
reordering of the different terms, and adding up, we 
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FIG. 2: Solid tilted and horizontal lines: Sketch of the Y- 
shaped SI— Q spectrum of p K . Dashed and dash-dotted curves: 
X-shaped SI — Q spectrum [Eq. for seed frequencies 

(circles) ±Sl s = ±0.63 (kr = 0.41) and for (squares) ±Sl s = 
±2.01 (kr = 2.23). 



readily obtain d^a = — «/ e d e a — (kr/SI s ) d T a at any dis- 
tance £, that is, Eq. ([5]). 

From the general solution of ([5]), the field at any 
propagation distance is found to be of the form a = 
F[p; ee~ K '£; r — (kr/SI s )£\ exp(ia£). In the asymptotic 
limit £ — * oo, ee~ KI ^ — > oo since kj < for an unstable 
perturbation, and the field must adopt the generic form 

a = G[p;T- tap(ia0 . (10) 

If as suggested by the observations, the nonlinear dynam- 
ics relaxes into a linear regime at large £, the field a in 
(JTUJ) can be identified with an X wave, or a diffraction- 
and dispersion-free solution to the linear Schrodinger 
equation in a medium with normal dispersion, [14| char- 
acterized by a inverse group velocity (or group velocity 
mismatch) Kr/SI s with respect to r = (i.e., with respect 
to a linear plane pulse at the Townes frequency). The X 
wave PH|) is equivalently described as a superposition of 
Bessel beams Jo(Qp) exp[i(Sl 2 /2 - Q 2 /2)£\ exp(-iOr) of 
different frequencies SI, and transversal wave numbers Q 
verifying the hyperbolic, or X-shaped, dispersion relation 

m 

Q= ^2[-a-(K R /Sl s )Sl + Sl 2 /2] (11) 

such that the axial wave number SI 2 /2 — Q 2 /2 — a + 
(k r /SI s )SI is a linear function of the frequency SI. 

Figure [2] shows two examples (dashed and dot-dashed 
curves) of the X-shaped SI— Q spectra to which perturbed 
Townes beams are expected to evolve when the Y pertur- 
bations are seeded at two particular couples ±Sl s of fre- 
quencies (circles and squares). The X spectrum presents 
always a frequency-gap, and does not keep any simple 
relation of symmetry with the seeded or Townes frequen- 
cies. The left branch of the X spectrum passes always 
close to the origin and intersects the tilted Y branch in 
the vicinity of —Sl s . The cut-off frequency of the right 
branch is always higher than +Sl s , approaching one to 
another as Sl s increases. 

Further algebra shows that the analysis leading to the 
X spectrum (|11[) would remain valid for Y perturbations 
at multiple couples of frequencies ±Sl s if the quotient 



hlr/SI s were constant at different perturbation frequen- 
cies Sl s . Though the value of k r /SI s obtained from the 
linear instability analysis is not constant, it presents a 
stationary point about the frequencies with maximum 
gain [see Fig. [Tja), dotted curve], by which it is reason- 
able to expect approximate X-like spectra also for multi- 
ple Y perturbations, as seen in Section HVl 



B. Numerical analysis 

We have tested the above predictions by numerical in- 
tegration of the NLSE |T]) by means of a standard finite- 
difference method, taking the perturbed Townes profile 

o(p,r,e = 0) = oo (p) + e[u(p)e-^ T + v*(p)e m « T ] 

= a Q (p)+p K (p,T,Z = 0) (12) 

as the initial condition. The initial spatiotemporal spec- 
trum a(Q, SI, 0) of a(p, r, 0) presents then only the three 
discrete frequencies SI = 0,±SI S corresponding to the 
Townes frequency and the Y perturbation p K . Fig- 
ure [3] shows [(a), solid curve] a(Q,Sl = 0,0) cx a (Q), 
or Hankel transform of the Townes profile, [(b), solid 
curve] a(Q,Sl = —Sl s ,0) cx ev(Q), and [(c), solid curve] 
a(Q,Sl = +SI S ,0) cx eu(Q), for particular choices of Sl s 
and e in the perturbation (| 1 2[) (see caption). 

On propagation the power in the Townes frequency 
SI = decreases [Fig. [UJa), dashed curves] and the per- 
turbation grows [(b) and (c), dashed curves]. Figure[3fd) 
displays the power redistribution among the different fre- 
quencies during the propagation. While the power loss 
in the Townes frequency SI = remains small (closed 
squares), the power in frequencies |0| > (closed cir- 
cles) grows exponentially with the gain — 2k/ (dashed 
straight line) as expected from the linear instability anal- 
ysis. The faster growth of the power in > observed 
at longer distances (£ ~ 8 and beyond) is seen to be 
due to the generation of additional couples of frequen- 
cies SI — ±NSl s , N = 2, 3, . . . with an exponential gain 
significantly higher than —2k\ (gray circles). These new 
frequencies extend over a spectral bandwidth that ex- 
ceeds by far the limit expected from the linear instability 
analysis (O s ~ 1.5). The spectra of radial frequencies 
Q of these new frequencies [see Figs. |U[a) and (b) for 
the particular choice SI — ±5Sl s ] feature a maximum at 
Q = 0, which represents a significant enhancement of 
the axial emission. More importantly, additional max- 
ima are observed to emerge at a values of Q that fits 
well, for each new frequency SI, to the dispersion relation 
(fTTj) of the X wave seeded at Sl s [vertical lines in Figs. 
SJa) and (b)]. Surprisingly, the X wave (fTTj) is seen to 
emerge well-before the relaxation of the dynamics with a 
strong exponential gain. 

Another unexpected feature is that the fast exponen- 
tial growth is accompanied by a slight drift of the X wave 
spectrum. The weak maxima of the Hankel transforms in 
Figs. BJa) and (b) fitting to (fTTj) with Sl s (vertical lines) 
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FIG. 3: Solid curves: Modulus of the Hankel transforms of 
(a) the Townes profile, and of the Y perturbation (b) at — £l s 
and at (c) +Q, S - The seed e is such that \a(p, T, 0)| differs 
from ao(p) by 0.1 as much, and the perturbation frequency 
is n s = 0.63. Dashed curves: (a) Depletion of the Townes 
Hankel transform, and growth of the Y perturbation at (b) 
— Q s and at (c) +£l s upon propagation, (d) In logarithmic 
scale, total power in the Townes frequency f2 = (squares), 
in the frequencies |fi| > (circles), and in the newly gener- 
ated frequencies > fi s (gray circles) versus propagation 
distance £. The dashed straight line represents the exponen- 
tial gain oc e" 2K/5 , with -2ki = 0.248 at Q B = 0.63. 



become gradually shifted to the sharper and much higher 
maxima in Figs. HJc) and (d). The shifted maxima fit 
also to the dispersion relation pip of an X wave (vertical 
lines) but with a perturbation frequency f2 s higher than 
that in the initial condition. The origin of the drift can be 
understood from Fig. [21 where we show X wave spectra 
(llip at two perturbation frequencies ±f2 s , corresponding 
to the firstly generated X spectrum (dashed curve) and 
the final, shifted X spectrum (dot-dashed curve) in the 
numerical simulation of Figs. [3] and 2] For the Y pertur- 
bation at ±f2 s in the initial condition (circles), the lower 
newly generated frequencies, say ±2f2 s , are characterized 
by a maximum at Q ~ v2|20 s | (enhanced conical emis- 
sion) for — 257 s , where the left X-branch and Y-branch 
are nearly indistinguishable, and by a single maximum 
at Q — (enhanced axial emission) for +2fi s . No addi- 
tional maximum appears at Q ^= since the frequency 
+2fl s is lower than the cut-off frequency of the right X- 
branch. The couple of new frequencies ±2f2 s then con- 
stitutes a secondary Y-perturbation to the Townes beam 
with slightly different quotient Kfl(2f2 s )/2S1 S , that gen- 
erates its own, shifted X wave. This process may be re- 
peated several times during propagation, giving the im- 
pression of a drifting X wave spectrum (jlip with some 
average Q s over all Y perturbations. 

Relaxation of the nonlinear dynamics occurs when the 
Townes beam has lost a significant part of its power 
(more than 50% in the example of Figs. [3J and 0|, 
and experiences therefore diffraction. As seen in Fig 
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FIG. 4: (a) |o(0 = -5Q S , Q, £ = 8)| (b) \a(Q = +5fi s ,Q,£ = 
8)|, (c) \a(n = -5tts,Q,£ = 16) | and (d) \a(Q = +5ft 3 ,Q,£ = 
16) | . The vertical dashed lines are placed at values of Q 
given by the X spectrum (111[) with perturbation frequency 
Q, s = 0.63 in (a) and (b), and with perturbation frequency 
Q s = 2.01 in (c) and (d). (e) Contour plot (10 decades in 
logarithmic scale) of the modulus of the spatiotemporal spec- 
trum \a(Q,Q,£ — 16) |. Dashed white curve: X spectrum 
(|lip with fi s = 0.63. Dashed black curve: X spectrum with 
Q s = 2.01. 



[3jd), interchange of power among the different frequen- 
cies ceases. Figure UJe) shows the final, asymptotic X- 
shaped spatiotemporal spectrum fi— Q, obtained by as- 
sembling the Hankel transforms at all frequencies ±NCl s 
(vertical strips) with significant power. The most ampli- 
fied spatiotemporal frequencies lie about either Q = or 
the X spectrum (jlip with the shifted £l s (black dashed 
curve) from the firstly generated X spectrum |[TT|) at the 
seeded perturbation frequency il s (white dashed curve). 
Extensive numerical simulations shows only quantitative 
differences in the reported dynamics for different values 
of the initial perturbation frequency Q s and its ampli- 
tude e <^ 1; no qualitative changes in the dynamics are 
observed. 



C. Four- wave mixing interpretation 

A simple picture of the above dynamics can be given in 
terms of the parametric amplification of new spatiotem- 
poral frequencies in a four-wave mixing interaction driven 
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by the Kerr nonlinearity, as suggested by the exponential 
growth of the new frequencies. 

The Townes beam is characterized by a frequency 
fl = and an axial wave number kx — a due to self-phase 
modulation [as in the NLSE (jTJ) , frequencies, wave num- 
bers and inverse velocities are referred to those of a linear 
plane pulse at the Townes frequency] . The Y wave p K is 
constituted by the two frequencies ±f2 s with wave num- 
bers k + = Q 2 /2 + 2a for +Q S (where the first term arises 
from normal group velocity dispersion and the second 
from cross-phase modulation), and fc_ = fig/2 — Q 2 /2 = 
—ft 2 /2 for — fl s (where the second term arises from tilt 
and the relation Q = v2| fi s | for Y waves has been used) . 
As a wave object, the Y wave itself represents an infinite 
train of pulses of carrier frequency [fi s + (— fi s )]/2 = 0, 
carrier wave number ky = (fc+ + fc_)/2 = a, beating fre- 
quency [fi s — (— ri s )]/2 = fi s , and beating wave number 
(k + — k-)/2 = Vl 2 /2 + a. Accordingly, the beatings of 
the Y wave p K propagate at a inverse group velocity (or 
group velocity mismatch) +fi s /2 + a/il s . 

Once the Y wave has grown up to a significant am- 
plitude, the Townes beam and the Y wave can act as 
pump waves, and amplify new, idler and signal waves at 
frequencies ±fi. Note first that a configuration in which 
the idler and signal frequencies ±fi form a new Y wave 
could be efficiently amplified since all Y waves of differ- 
ent frequencies ±fi are phase-matched with the pump 
Townes and Y wave. However, the beatings of the new 
Y wave propagate with different inverse group velocity 
+fi/2 + a/fi from those of the pump Y wave at ±fi s . 

More efficient amplification is expected for idler and 
signal frequencies that are not only phase-matched but 
also group-matched with the pump waves. For idler wave 
at +fi with k\ = fi 2 /2 — Q 2 /2, signal wave at — fi with 
/c2 = fi 2 /2 — Q?,/2, producing beatings at inverse group 
velocity (— Q\+Q 2 ,) /2Q,, the condition of phase-matching 
(k T + k Y = fci + k 2 ) reads as 2a = fi 2 - (Q 2 + Q|)/2, 
and the condition of group-matching as fi s /2 + a/£l s = 
(—Qi + Ql)/-^- These two conditions lead to 



Q1.2 = 



-a 



n + - 



(13) 



where Q\ is obtained for positive fi and Q2 for negative 
fi. Equation (p~3|) is the same as the X-shaped spectrum 
([TTj) . provided that ~ fl 2 /2+a, which in fact has been 
shown in [6| to approximate kr where the gain —2kj is 
significant [see dot-dashed curve in Fig. Ufa)]. 

The X-shaped spectrum is then explained as the re- 
sult of the preferential amplification of phase- and group- 
matched spatiotemporal frequencies by the Y-perturbed 
Townes beam. In particular, the group velocity of the 
X wave is the same as the group velocity of the beat- 
ings of the Y wave. Due to the nonsymmetric frequency- 
gap in the X spectrum, however, the optimum condi- 
tions of phase- and group-matching cannot be attained 
for ±il close enough to 17 = 0, in which case a new phase- 
matched Y wave is amplified. 
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FIG. 5: Dashed curve: weakly perturbed Townes beam \a(p = 
0,t, £ = 0)j as given by Eq. (I15[) . The amplitudes of the 
two Y pulses are chosen such that e KjS = e_ K * jS and to yield 
|p«, s (0, 0,0)| = |p_ K *, s (Q, 0,0)| = 0.0025 at all frequencies Q s 
in the interval [0, 1]. For this choice, the two weak Y pulses 
overlap initially. Solid curves: amplitude \a(p = 0, t, £)| at 
increasing propagation distances £. 



IV. NONLINEAR SPATIOTEMPORAL 
INSTABILITY: DOUBLE Y PERTURBATIONS 
AT MANY FREQUENCIES 

The preceding analysis allows us to understand the 
more realistic situation in which the self-guided beam 
experiences perturbations at a continuous of frequencies 
±J7 S about £1 = 0. To this purpose we have simulated nu- 
merically the evolution under the NLSE (fTJ) of an initial 
condition of the type 



o(P)T)C=0) = a a(j>) + 



(14) 



+ / dn s e K , s [u s (p)e- in ° T +v*(p)e m ° T ] 

JQ 



+ / dn s e^ K . iS [<(p)e- insT +u s (p)e^ r ] 
Jo 



or, equivalently, 



a(p,T,£=0) = a (p) + / dn s p K . s ( P ,T,^0) (15) 
Jo 



+ / dfl s p- K *. s (p, t, £ = 0) , 
Jo 

that describes the Townes beam perturbed by the two 
unstable Y waves excited at a many frequencies ±f2 s 
with small amplitudes e K , s and e- K * jS . Each integral in 
(fTS")) will be referred to as a weak Y pulse, since it repre- 
sents a perturbation of finite duration with the Y-shaped 
spatiotemporal spectrum of Fig. QJb) (and its reflection 
about Q s = 0) and small amplitudes determined by e KiS 
(and e_ K * jS ) at each frequency ±fi s . An example of this 
type of perturbation, with frequencies Q s filling the in- 
terval [0, 1] is shown in space-time domain in Figs. [5] 
(dashed curve) at the Townes beam center p = 0, and in 
Fig. E](top) for p > (see captions for numerical values). 



^■3ts lacal time r 




lacal time t 

FIG. 6: Top: weakly perturbed Townes beam \a(p, t,£ = 0)| 
as given by Eq. (|15[1 . Numerical values are the same as in 
Fig. [5] Bottom: modulus of the amplitude \a(p, t, £ = 13.75) 
of the propagated field for the initial condition on the top of 
the figure. 



Numerical integration of the NLSE (fTJ) with a only one 
weak Y pulse perturbation (e_ K * s = or e K s — 0) shows 
that it develops in the nonlinear regime of instability into 
a single, steepening, subluminal pulse (when e_ K *. s = 0) 
or superluminal pulse (when e KtS = 0) of amplitude ex- 
ceeding the Townes amplitude. If instead the two weak Y 
pulses are excited (e_ K * jS ^ and e KiS ^ 0), two steepen- 
ing, subluminal and superluminal pulses grow, as seen in 
Figs. O (solid curves) and [5] (bottom). These two giant Y 
pulses can then grow independently and result from the 
exponential growth of the two weak Y perturbations in 
Eqs. ([Fiji or (fT5)) beyond the linear regime of instability. 
Thus, with little loss of generality, we have limited the 
numerical simulations to the case e K . s = e_ K * iS leading 
to symmetric spectra about O = and symmetric pulse 
shapes about r = 0. 



Subluminality and superluminality, and self-steepening 
on opposite fronts during propagation, are distinctive fea- 
tures of the Y pulses, that originate from their reflected 
Y-shaped spatiotemporal spectra. These features are il- 
lustrated in Fig. [7] by means of a simplified model of 

Y pulses. Contrarily to a standard Gaussian pulse, the 
positive and negative frequencies in the spectrum of each 

Y pulse experience group velocity dispersions of opposite 
signs. For the Y pulse e.g., the frequency +f2 s prop- 
agates along the ^-direction subject to normal dispersion 
(k + = Qg/2+2a), and the frequency — ft s does subject to 



FIG. 7: Simplified, one-dimensional model of propagating Y 
pulses to evidence the opposite group velocities and steepen- 
ing in opposite fronts. A pulse at £ = (dashed curves) with 
a Gaussian spectrum [exp(— Q 2 /AQ 2 ), with Af2 = 3] is made 
to propagate with anomalous (normal) group velocity disper- 
sion for positive (negative) frequencies in (a), which models 
P- K *, s , and vice versa in (b) for p re , s . Solid curves: £ = 0.4, 
Dotted curves: £ = 0.9. 



anomalous dispersion (fc_ = —Q^/2), this being a con- 
sequence of its tilted propagation direction. Thus, dif- 
ferently from a Gaussian pulse, opposite portions of the 
spectrum travel at identical inverse group velocities along 
the ^-direction (dk+/dQ s \n s = dk- /dQ s \-n s = +f2 s for 
p K)S ), which prevents from any pulse broadening associ- 
ated to pulse chirping. In particular, if ±O s is a suit- 
ably defined mean frequency of the positive or negative 
half-spectra of the Y pulse p K . Sl then its group velocity 
as a whole is determined by +Cl s , i.e., p KjS is sublumi- 
nal. Similar considerations for the Y pulse p_ K * jS lead 
to superluminality with inverse group velocity — Q s . The 
opposite group velocities of the two Y pulses can be seen 
in Fig. where the two halves of a Gaussian-distributed 
spectrum are forced to propagate with opposite group 
velocity dispersions. 

Though Y pulses remain unchirped during propaga- 
tion, they experience broadening and distortion (see Fig. 
UJ by a different mechanism from that of standard pulses. 
Broadening of Y pulses is caused by a chirp in the beat- 
ings between each couple of frequencies ±J7 S due to a dis- 
persion in the beating group velocity, i.e., to a frequency- 
dependent group velocity f2 s /2 + a/VL s for the beatings 
in the Y pulse p KlS , and — f2 s /2 — a/Cl s for beatings in 
the Y pulse p- K * tS . Moreover, assuming a well-behaved, 
compact spectrum about O = with mean positive and 
negative frequencies ±fi a , the frequencies ±217 s define 
roughly the extreme frequencies with significant ampli- 
tude in the spectrum, and therefore the fastest beat- 
ings in time. They propagate at inverse group velocities 
2fi_ s /2 + a/2n s ~ +Q S for p K , s , and ~2Q s /2 - a/2fl s ~ 
—fl s for p_ K * iS (where the approximation holds for rel- 
atively narrow spectra, fl s ^> y/ot/2 — 0.32). In other 
words, the fastest beatings travel at group velocities ap- 
proximately equal to the group velocities of the Y pulses, 
forming thus the steep fronts about the center of the Y 
pulses (see Fig. [7J, whose rise times are seen to be di- 
rectly linked to a frequency ~ 20 s . 
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£ =15 




FIG. 9: (a) In logarithmic scale (9 decades plotted), contour plot of modulus of the spatiotemporal spectrum a(Q,il,£ = 6.5) 
for the initial condition in Figs. [5]and|5l where the two Y pulses were seeded at frequencies fl 3 in [0, 1]. (b) Fitting double-X 
spectrum of Eqs. (|16[) an d (|17|l with £l s = f2 a (£ = 6.5) ~ fi s (£ = 0) = 0.5 (broadening is very small), (c) In logarithmic scale 
(9 decades plotted), contour plot of the modulus of the amplitude of the spatiotemporal spectrum a(Q, fi, £ = 15) for the same 
initial condition, (d) Fitting double-X spectrum of Eqs. (fl6)) and (fT7|) with fl s = 2f2 s (f = 15) = 5.97. 




FIG. 8: In logarithmic scale (9 decades plotted), modulus of 
the spatiotemporal spectrum a(Q,fl,£ — 16.25) of the prop- 
agated field a(p, r, £ = 16.25) for the initial condition in Figs. 
[S]and[51 Solid and dashed straight lines: double Y spectrum. 
Solid and dashed curves: double X spectrum. 



The actual shocks in the simulation of Figs. [5] and 
[5] show an unbounded increase in the slopes, which can 
be explained as the combined effect of Y pulse steepen- 
ing and of increase of the extreme frequencies, that is, of 
spectral broadening in the form of new, phase-matched Y 
waves at frequencies beyond those in the initial condition, 
as explained in Sec. IIIII for single Y wave perturbations. 
Figure [5] displays the hot portion of the spatiotemporal 
spectrum of the propagated field for the initial condition 
in Fig. [6] (top), and evidences a significant Y-shaped 
spectral broadening. The double Y spectrum, character- 
ized by maxima at the straight lines Q ~ \/2|f2 s | and 
Q = 0, is broadened in the propagated field well-beyond 
the interval [—1, 1] in the initial condition. The fine ver- 
tical structure observed in this spectrum originates from 
interference between the two splitting Y pulses, and the 
broad maxima above the Y-shaped spectra are portions 
of two X-shaped spectra, as described below. 



The nonlinear interaction between each Y pulse and 
the Townes beam results in the amplification of the two 



X waves of opposite group velocity mismatches 



y/2[-a- (kr/O s )0 + 2 /2], (16) 



8 



Q p _ H , t , = y/2[-a + (K R /n s )n + Q 2 /2] , (17) 

as seen in Fig. [5J or in Fig. [5]spanning a much wider spa- 
tiotemporal frequency window. The process of amplifica- 
tion of each X wave from each Y pulse follows the same 
trend as for single- frequency Y wave perturbations, as ex- 
pected from the slow variation of the inverse group veloc- 
ities, kr/£I s ~ Q s /2+a/Q s or —kr/£I s ~ — Q s /2— a/Q s , 
across the spectrum of each Y pulse. Amplification of the 
two X waves are independent processes one from another 
due to the non-overlapping spatiotemporal regions of in- 
teraction. In fact, the two X-shaped spectra in Eqs. (fl"6|) 
and (fl~7|) are seen to develop only when the two Y pulses 
have split-off. Figure [9^a) shows an example of the in- 
cipient, double- X-shaped spectrum that emerges immedi- 
ately after Y pulse splitting (see also Fig. [5] at £ = 6.5). 
The location of the maxima of this incipient double-X 
spectrum is seen to fit accurately to Eqs. (fl"6| and (fTTjl . 
as seen in Fig. OJb), with a perturbation frequency fl s 
equal to the mean frequency f2 s (£) of the half-spectrum, 
defined by 

o / oo dn/ oo dQQ|a(Q,n,Q|n 

" slU " f~<KlJ~dQQ\a(Q,n,0\ ■ 

With the development of the shock fronts due to the gen- 
eration of new phase-matched Y waves, the two growing 
X spectra start to drift apart, as described for single- 
frequency Y wave perturbations. Figure |^c) shows a 
shifted double-X spectrum at a distance where the shocks 
are well-developed, and Fig. [9jd) its fitting to Eqs. (fT6|) 
and (|17| . Noticeably, the drifting X waves are no longer 
described by Eqs. (fT6|) and (fT7|) with the "instantaneous" 
mean frequency f2 s (£) of the broadening spectrum, as 
could be expected, but with il s close to 2fi s (£), i.e., twice 
the instantaneous mean frequency. 

This fact can be interpreted to be the result of the pref- 
erential amplification of X waves that are highly coupled 
to the Y pulses and the self-guided beam. Each X wave 
is 1) phase- matched with each Y pulse and the Townes 
beam, as explained in Sec. II VC . 2) group-matched with 
the fastest beatings at 2f2 s in each Y pulse forming its 
shock, as also explained in Sec. IIVC and in relation to 
the simple model of Y pulse. In addition, each X wave 
is 3) group-matched with each Y pulse, since the group 
velocity of a Y pulse coincides with the group velocity of 
its fastest beatings, as seen in the model of Y pulses. 

V. CONCLUSIONS 

Though an accurate description of axial spectral 
broadening, off-axis spectral broadening or conical emis- 



sion, pulse splitting, shock formation, and development 
double- X-shaped spatiotemporal spectra with frequency- 
gap s in the propagation of a light filament, 0, la H, El 
llll.fI2l.fT3T] requires the inclusion of all relevant, linear and 
nonlinear material effects, they have been found to exist 
with essentially the same characteristics in the spatiotem- 
poral dynamics of a self-guided light beam in a Kerr 
medium with normal group velocity dispersion. The key 
concept unifying all these phenomena is that of unstable 
Y-shaped wave modes, that can grow from infinitesimal 
perturbations to the self-guided beam, and are seen to de- 
velop into two giant, spectrally broadened, splitting and 
shocking pulses, amplifying two super-broadened, phase- 
and group-matched X waves. 

In particular, the shock fronts in the trailing edge of 
the subluminal pulse and in the leading edge of the super- 
luminal pulse, have been discovered recently in [l3T ] , and 
have been described as shocking-X waves arising from 
the subluminal (or supcrluminal) propagation and spec- 
tral broadening. These shocks appear here as the result 
of linear Y pulse steepening and spectral broadening. 

Inclusion of higher-order dispersion effects is straight- 
forward and would lead to the so-called "fish-like" wave 
modes, or distorted X waves, in place of the hyperbolic X 
waves, as observed in experiments and numerical simu- 
lations in conditions of non-negligible high-order disper- 
sion, pi]. 

In comparison with previous theories of nonlinear X 
wave (or more generally, fish-like wave) formation, the 
splitting pulses and their group velocities are not taken 
here from numerical simulations of the full Maxwell equa- 
tions including all material effects, [ll], [l^] or are not 
regarded as a source terms with prescribed group ve- 
locities in an essentially linearized Schrodinger equation 
[l6| . Splitting with a certain group velocity mismatch ap- 
pears here as an intrinsic feature of the unstable Y modes 
that can grow from (non-splitting) perturbations at noise 
level. These alternative theories have been seen to admit 
also a simple interpretation in terms of a three- wave mix- 
ing interaction (rather than four- wave mixing), (lH[l3| in 
which an X wave results from the scattering of an incident 
wave by a split-off (nonlinear polarization) pulse, which 
predicts a linear growth of the X wave, [16| against the 
exponential and growth obtained from the present the- 
ory. Further experiments could elucidate these points. 
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